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We consider a bouncing Universe model which explains the flatness of the primordial 
scalar spectrum via complex scalar field that rolls down its negative quartic potential and 
dominates in the Universe. We show that in this model, there exists a rapid contraction 
regime of classical evolution. We calculate the power spectrum of tensor modes in this 
scenario. We find that it is blue and its amplitude is typically small, leading to mild 
constraints on the parameters of the model. 



1 Introduction 

Inflation is the best studied and very plausible scenario for the cosmological epoch preceding 
the hot Big Bang era. Nevertheless, it is of interest to consider alternative scenarios. One 
of them is Ekpyrosis [U [2] , understood in a broad sense as a period of slow contraction with 
equation of state p > p- Its particular version, pseudo-conformal Universe [3l H] makes use of 
a nearly conformal evolution to generate scalar perturbations with nearly flat power spectrum, 
thus building upon the idea that the flat spectrum may be due to conformal invariance [SI El [7]. 
The adiabatic perturbations generated by conformal mechanisms have interesting non-linear 
properties, such as non-Gaussianity and statistical anisotropy [8]. 

The concrete example of the pseudo-conformal Universe [3] is based on the Einstein gravity 
interacting minimally with a complex scalar fleld with negative quartic potential. The action 
is (signature ( — h ++)) 



M2 

S = — I d'^Xy/^R + / d'^Xy/^ 

lOTT 



The contracting epoch is assumed to begin with the asymptotically Minkowski space-time and 
= 0. Since the potential is negative, this state is unstable, and under the assumption of the 
homogeneous background evolution, the fleld rolls down as follows. 
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At this stage perturbations of the phase Arg are generated, which, to the hnear order, 
automatically have flat power spectrum [6]. These entropy perturbations are assumed to get 
reprocessed into adiabatic ones by one or another mechanism, see, e.g., Refs. O [TOl [11]. It 
has been pointed out [3] that tensor perturbations are generated as well, and that they have 
blue power spectrum. The latter is similar to the tensor spectra generated in other ekpyrotic 
models fT^. 

In this paper our focus is on the late-time contraction in this model. The classical behavior 
turns out to be rather interesting: the initial slow contraction with the effective equation of state 
p ^ p is followed by the rapid contraction regime in which the equation of state asymptotes to 
p = p. As in all ekpyrosis-like models, one assumes that the latter regime ends up in a bounce 
and defrosting [13] whose mechanism we do not specify. 

We also consider the generation of tensor modes, both at slow and rapid contraction stages. 
We find that their power spectrum is always blue, and the amplitude is small enough, so that 
no strong constraints on the parameters of the model can be placed. Thus, the model is safe 
from this viewpoint. 

This paper is organized as follows. In Section 2 we consider the classical contraction. In 
Section 3 we study the tensor perturbations. We conclude in Section 4. 



2 Classical solution 

The homogeneous contracting Universe is governed by the equation of motion for the field 0, 

+ 3H<i) = X(f)^ , (2) 

and the Friedmann equation 

r2 Svr /I - 2 A, ,,4 



At early times, t — )■ — cxd, one neglects the Hubble friction term in the first equation and obtains 
the solution ([T]), where we take real. To this order, the energy density of the field van- 
ishes. The Hubble parameter and then the scale factor can be obtained from the Raychaudhuri 
equation 

2d [dV Svr / 0^ A0^' 



One finds 

Stt 1 

3AM2(-t)3- (4) 

This corresponds to the equation of state 

p 



—7- oo 

P 



From eq.(|4]) we get the scale factor 



where we keep the overall normalization factor fl** as a free parameter. So, the Universe 
contracts slowly at early times. 
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We are mainly interested in the solution at late times. For large H we neglect the scalar 
potential and solve the simplified equation of motion and Friedmann equation: 



(j) + 3H(j) = 

Upon substituting H into the first equation we get 

•■ _ V^127r -2 



(6) 



pi 



The solution is 



pi 



ln(/i(t, -t)) 



(7) 



where and yU are to be determined bymatching this solution to eq. ([T|. We can indeed 
neglect the scalar potential in eqs. ([2), (3|, provided that Xcp"^ ^ 0^ and X(p^ <^ 0. Using eq. 
([T]) we see that these conditions are satisfied at late times, such that 



t < 



127r 



A Mpz|ln(/x(t, -t))|3/2 



Note that for small A, there exists a time interval where the evolution is still classical, (t* — t) ^ 
Mp^\ while the inequality (8) is satisfied. 

It follows from eqs. ^ and ([T]) that at late times, the Hubble parameter and scale factor 

are 

o = <(i. -^)"^ (10) 

where a'^ is yet undetermined parameter. They correspond to the equation of state 

P = P ■ 

Let us match the two regimes and express the coefficients t^:,a'^ and /i in terms of a=„^,. For 
estimates, we require that 0, a and their time derivatives are equal for the solutions ([T]), (jsj 
and ([T]), ([To]) at some time tjv/- We write 



0: 



4> 



H : 



V2 



^(-tA/) 

V2 



'X{-tM? 



in 



pi 



Vl2^ 
1 



In -Im)) 



pi 



VT27r — tM 



,1/3 



3AM2( 



-tM) 



3(t* — tM) 



(11) 

(12) 

(13) 
(14) 



As a cross check, the potential and Hubble friction terms in eq. ^ are of the same order at 
time tM- A0^ ~ 3H(p. Likewise, the kinetic and potential terms in the Friedmann equation ([3) 
are of the same order at that time, 0^ ~ A0^. 
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By solving eqs. ( 11 ) - ( 14 ) , we obtain the time of the transition from slow to fast contraction, 



/8vr 1 
V 3 VXMpi 

Of course, if the bounce happens at t < tM, the fast contraction regime does not occur at all. 
Assuming that this is not the case, we find the parameters of the fast contraction stage 

/i ~ VXMpi (16) 



/327r 1 



An example of numerical solution is shown in Fig. [T| One observes that the two regimes ([T]) 
and ([T]) indeed nicely match at intermediate time. 

In what follows we use the solution in the second regime, written in terms of conformal time 
(in the first regime conformal time equals t/a**): 

a = a,ir],-7]y/\ (19) 



Making use of eqs. (17), (18) and requiring that r](tM) = tM/O'**-, we get 



al ^ alMpiVX , (20) 



(21) 

To end up this Section, we express the parameters of our solutions in terms of physical 
quantities. We consider the simpliest possibility that the contracting stage terminates abruptly 
by the bounce and instantaneous defrosting. We therefore assume that both the scale factor 
and energy density (and hence \H\) coincide at the end of the contracting stage and at the 
beginning of the hot epoch. At the beginning of the hot epoch just after the bounce we have 



7^2 

A 
hi 



Hb = T7^ (22) 
M*, ^ ' 



1/3 

(23) 



where is the defrosting temperature, M*i = M*i{Td) = i ee^q^d refers to any moment 

of time at the radiation dominated epoch. On the other hand, just before the bounce the 
Hubble parameter and scale factor are given by 

We thus obtain 

1 m;, M;,gyf 

r/, - r/6 = — — f- = — , (25) 

2a6 TJ 2aTTTdgl!^ 



4 







Figure 1: An example of numerical solution to eqs. (|3| (solid line). Short dashed line is 
^/2/{-^t). Long dashed line is -4= ln(/i(t* - t)). 



1/2 



Finally, the scale factor at the beginning of contraction is 



(26) 



1/3 



(27) 



Let us now turn to the calculation of the power spectrum of tensor perturbations. 



3 Tensor perturbations 

In this Section we work in terms of the canonically normalized field 



'Ml 



3277 



(2^ 



where h is the shorthand for the transverse traceless component of the metric perturbation. 
The field equation has the following form in conformal time: 



+ 2-(^' + fcV = . 
a 

In terms of the variable x = this equation reads: 

„ a" 



^2 



X 



-x + = . 



(29) 



(30) 
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In what follows we consider two regimes separately: early times, when the scale factor a is 
varying slowly; and late times, when it evolves as a oc (t* —tY^^. The former regime is relevant 
for low momentum modes which exit the horizon early, while the latter is important for high 
momentum modes. 

We begin with the modes of low momenta. In the slow contraction regime, the scale factor 
behaves as follows: 



a = a* 



7]^ 



An 



The low momentum modes are those which exit the horizon, k/a ~ H, well before time tM- 
Making use of eqs. ^ and (15), we obtain for these modes 



Note that for these modes one has krjx <^ 1, where rjx is the conformal time of horizon exit. 



(31) 



In the slow contraction regime, eq. (|30|) has the following explicit form: 

. 



6c 

X" + ^X + k\ 



(32) 



We assume that the field is in its vacuum state at past infinity and hence has the standard 
asymptotics x^"^^ = g±«fcf? q^^^ large After crossing out the horizon, but before the onset of 

rapid contraction, the solution is 



X = Cia + 



dt 



7] 



V 



3c 
V 



(33) 



Now we match the two asymptotics and find x^~^^ ^iid x^~^ iii the superhorizon regime. To 
this end we solve eq. (32) perturbatively near horizon crossing. In the slow contraction regime 
we have ^ <^ 1, and the horizon crossing occurs when /c^ ~ ^. Hence, there are two small 
parameters around the horizon crossing: k'^t]'^, ^ <C 1, and the two linear independent solutions 
to eq. (32) can be found as series in these parameters: 

1 



X{i) 

X{2) = V 



7]^ 

3c 



(34) 
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Before the horizon crossing, when ^ <^ k'^f]'^ <^ 1, these solutions match to X{i) = cos{kri) and 
X(2) = sm{k7])/k. Since x^ 



/(+) gg^ 



X^^^ = X(i) + ikx(2) ■ 



After horizon crossing, when k rj <^ ^ <^ 1, this solution reads 



X 



(+) 



+ ik [t] 



3c 

7] 



(35) 



which coincides with eq. (33) with Ci 
horizon exit is 



( + ) 




i/ca**. So, the low momentum mode after 



(36) 



where we inserted the standard normahzation factor and made use of the fact that fc?7x ^ 1. 
We cross check this calculation in Appendix. 

The result (36) is valid also at late times, when the Universe contracts rapidly. Making use 
of eqs. ([5]) and (19), we calculate the integral in eq. (36) and find that in the rapid contraction 
regime, the low momentum modes behave as follows: 



(+) 



(27r)3/V2A; \ a** 



6c 



A;2 



1/4 



1 tu 



2c 



2c fk 



M) 



2 \ 1/4 



6c 



In 



V* 



It follows from eqs. (15), (20) and the condition of smallness of momentum (31) that the main 
contribution comes from the first, constant term: 



{+) 



1 



1/3 



(27r)3/V2A; \ 
With the standatd definition of the power spectrum. 



{ip{x,t)(p{x,t)) 



Ank^ 



we get 



4^ 



2/3 



(37) 



So, the power spectrum is blue, in agreement with Ref. [3j. According to eq. (28), the power 
spectrum of metric perturbations is 



_ 9 327r _ 16P 



2/3 



where the factor 2 accounts for two polarizations. 

Let us now turn to high momentum modes, such that k ^ \/~\Mpia^^, cf. eq. (31). These 
modes exit the horizon at the rapid contraction stage, so the field equation is 



X + 



1 



X 



4 [r]^ - riY 

Its solution that tends to e*^'' at large negative t] reads 



+ fc\ = . 



X 



(+) 



nk 



iv* -V)H'i^ [{V* -ri)k] ■ 



After horizon crossing, this mode has constant part and logarithmically growing one. The 
dependence on k is only logarithmic: 



{+) 



1 



(27r)3/V2A; a 
Hence, the power spectrum is again blue. 



^(+) 7r + 2zln(^i^) +2^7 



(3^ 



Kik) 



1 + 



7 + In 



k{ri^ - f]) 



(39) 
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This power spectrum refers to time before the bounce. To obtain the power spectrum after 



the bounce, we match the solution (38), vahd before the bounce, to the general superhorizon 



solution after the bounce, i.e., at radiation domination, when a{r]) = a{r] — rj). Using the 



instantaneous bounce approximation, cf. eqs. (24) - (25), we get 



Vb-V 



ajflb) 



2(^* - Vb) 



The superhorizon solution to eq. (29) at radiation domination is 

B 



W = A + 



1] — 1] 



By matching this expression and its first derivative to eq. ( |38[ ) at r/ = rjb, we get at late times 

TT + 2i(7 - 2) 



4^20, 7r2 



So, the primordial power spectrum at radiation domination is 



I.e., 



1 + 



8k^ 



1 + 



7 - 2 + In 



7 - 2 + In 



k{v* - Vb) 



(40) 



Let us now discuss possible effects of the tensor modes in the present and recent Universe. 
To this end, we again recall our assumption of the instantaneous bounce and defrosting. We 



make use of eq.(27) and write the present value of the critical momentum, eq.(31), as follows 

To 



= ^/XMpia 



'g*d 



where the subscript refers to the present value, and we set = 1. We see that this momentum 
is fairly high for not too small A. 

We first consider direct search for gravity waves. Barring the very small value of A, the 
relevant wavenumbers are below the critical value, k < kc- Hence, these modes exit the horizon 



at the slow contraction stage, and we use eq. (37) for their power spectrum. We again recall 



eq. ( 26 ) and find for the energy density 

dpGW 



dink 



2 ■ {2TxYeV^{k)ai = k 



>TiTTx^l^gll^gl',\l.m) 



2/3 



pi 



1/3 

9^k 



2i 



where is the scale factor at horizon re-entry at radiation domination, = ( 9*t_ 



(41) 



1/6 



The result (|41|) is hopelessly small: as an example, for k ~ 100 Hz we find that even with 

dVtcw 



T(i ~ Mpi and A ~ 1 



din A; 



5 X 10 



-24 



which is far below the projected sensitivity of gravity wave detectors. 
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We now turn to the BBN constraint 

PGwitNs) < ^r^5 • (42) 

We consider the case in which the rapid contraction regime does occur at the contracting stage 
(in the opposite case one obtains even weaker constraint on the parameters of the model). Since 
the power spectrum is blue, we make use of the result (40) and obtain 



J \aNs/ a^gk M , 



2 r)^2 rp4rp4: 

^pl 

where kmax is the highest momentum of modes that ever exit the horizon, kmax/^b ~ -fffe, or 
kmaxiv* ~ Vb) ~ 1- Thus, the constraint (42) translates into a weak bound on the defrosting 
temperature 

T. < 0.22 ■ 



1/3 1/6 



Thus, the pseudo-conformal Universe model is viable in a wide range of its parameters. 



4 Conclusion 

To summarize, we have solved the classical equations of motion and found the evolution of the 
Universe in the pseudo-conformal model. There are two regimes, one after another. The first 
is ekpyrotic (slow contraction) and the second is fairly rapid collapse. We then calculated the 
tensor spectra generated in the two regimes. We have seen that in both regimes there is no 
considerable production of tensor modes, hence this model is acceptable from the viewpoint of 
generating the cosmological perturbations. 
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Appendix 



To cross check the calculation leading to the result (35), we explicitly solve eq. (32) in the limit 
^ 1. We write 

X = e'''^{l + f{v,k)). 



Then eq. (32) takes the following form: 



6c , 



f" + 2^kf' + -(l + f) = 0. 



(43) 
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For c/rf ^ 1, the correction is small, / <^ 1, so that eq. (43) is 



The solution to this equation, which asymptotes / = as ?7 — )■ — oo, is 

V v' 

I 

—oo — oo 

Upon changing the order of integration we get 

— oo 

and after three integrations by parts we come to 

^2 ^ J rj' 

— oo 

This means that in the limit krj <^ 1 the solution is 

ikn A c 2ikc \ ^ c 3ikc 



7]"^ 7] J 7]"^ 7] 



This confirms our result (35). 
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